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Abstract 

We study the tunneling into the edge of the Luttinger liquid with periodically 

modulated concentration of electrons. It is shown that the modulation, by 

causing a gap in the spectrum of plasmons, leads to an additional anomaly in 

the density of states at frequency corresponding to the center of the gap. The 

shape of the anomaly depends strongly on the phase of the modulation. The 

sensitivity to the phase is related to the plasmon mode, localized at the edge, 

its frequency lying within the gap (analog of the Tamm state for an electron 

in a periodic potential) 

PACS Numbers: 73.20.Dx, 73.40.Gk 



1 



Typeset using REVIgX 



I. INTRODUCTION 

The explicit solution of the model of one- dimensional interacting electron gasS (Luttinger 
liquid) allowed to describe analytically the gap in the tunnel density of states for the region 
of energies close to the Fermi level. The physical origin of the gap is that, in order to add 
an electron to the interacting system, a shift, vacating an extra space, should be created. 
Then the gap reveals the smallness of the overlap of many-electron wave functions with- 
and without such a shift. The bosonization transformationEH! permits one to present the 
formation of the shift as a creation of plasmon modes in the system. The fact that the 
dispersion law of plasmons is linear in one dimension: uj = sk, results in a power-law 
behavior of the density of states v{u). For tunneling into the edge of a single-channel liquid 
the density of states has a formBi: 

s 

v{u) (XU) v f , (1) 

where Vp is the Fermi velocity. Matveev an Glazmani have traced the evolution of the 
exponent in ([[]) with increasing number of channels. 

Bosonization procedure can be generically extended to the the case of longitudinally 
inhomogeneous Luttinger liquid. If the spacial scale of inhomogeniety is much larger than 
the Fermi wave length, kp 1 , it can be considered as a "scattering potential" for the plasmon 
modes, so that a plasmon gets partially reflected after propagation through the region of 
inhomogeniety. If the boundaries of this region are abrupt compared to the wave length 
of a plasmon (but smooth on the scale of kp 1 ), the reflection coefficient, as a function of 
frequency, exhibits the interference oscillations. Such oscillations were first considered in 
Ref. |6|. In Ref. |7| it was demonstrated that the oscillations of the reflection coefficient lead to 
the oscillating structure in the density of states, superimposed on the power-law increase. 

In the present paper we consider another realization of inhomogeniety. Namely, we 
assume that the density of the Luttinger liquid is periodically modulated with a period 
a ^> kp 1 . We show that such a modulation results in additional singularity in the tunnel 
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density of states. The origin of the singularity is a gap in the plasmon spectrum, which 
opens due to the periodic modulation. In complete analogy with motion of an electron in a 
periodic potential, the position of the center of the gap, ujq, is determined by the condition 
k(uo) = 7r/a, so that uq = its /a. We assume the modulation to be weak. Then the width 
of the gap, A, is much smaller than u , the ratio A/ujq being proportional to the relative 
amplitude of the modulation. It is also proportional to the interaction strength, since it is 
obvious that without interaction the gap is absent. In the vicinity of ujq the modulation- 
induced correction to the tunnel density of states, 5u, is a function of uj = 2{oj — uo)/A. 
For \uj — ujq\ ^> A, 5u(uj) behaves as a power law with the same exponent as in ([[]). Since 
we study the tunneling into the edge, the actual shape of 5u(u) appears to depend strongly 
on the phase of the modulation. This is because for plasmons with frequencies close to u 
the behavior of the field near the edge is determined by this phase. The other reason is 
that, with periodic modulation, the presence of an edge leads to the formation of a localized 
plasmon mode with frequency within the gap. For an electron in a periodic potential this 
fact was established more than 60 years ago!. The frequency of the localized mode sweaps 
through the gap as the phase of modulation changes. We show that the localized plasmon 
provides a contribution to 5v{uj) comparable to that from the extended modes. 

The paper is organized as follows. In the next section the form of the Hamiltonian for an 
inhomogeneous Luttinger liquid is established. In Sect. Ill we derive a general formula for 
inhomogeneity-induced correction to the tunnel density of states. In Sect. IV the plasmon 
modes for the case of a weak sinusoidal modulation of the electron concentration are found. 
The results for 8v{u) are analyzed in Sect. V. Section VI concludes the paper. 

II. HAMILTONIAN 

In this section we will establish the form of the Hamiltonian for the Luttinger liquid with 
spatially varying concentration of electrons uq(x). We start from the classical equation of 
motion 
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where m is the electron mass; u(x,t) is the displacement; n(x,t) is the electron density: 
n(x,t) = rio(x) + nx(x,t), with n\ describing the temporal fluctuations; S(x,t) = £q(x) + 
£i(x, t) is the actual electric field consisting of static and fluctuating parts; P = -K 2 h 2 n 3 /3m 
is the hydrostatic pressure. The static electric field can be eliminated from Eq. (1) using 
the equilibrium condition 

dP 



e£ n = - 

ax 



(3) 

n=no(x) 



which gives £q{x) = — (ir 2 h 2 n /em)(dno/dx). Substituting this value into (1) and linearizing 
(1) with respect to u, £\ and ni, we obtain 

d 2 u TT 2 h 2 ( 2 dn\ dno\ 

mn — = -en o£l - — [n — + n^—j . (4) 

Let V(x) be the potential of electron-electron interactions, screened by the presence of 
a gate electrode. Then £\ can be presented as a field created by the density fluctuations 

n\(x, t) 

d r 

e£ x = — - / dx'Vix - x')n x (x', t). (5) 
dx J 

As usually we will assume that the screening radius is much smaller than the characteristic 
spatial scale of the fluctuations. Then Eq. (Q) reduces to e£\ = —V (dni(x,t)/dx), where 
V = JdxV(x). 

As a next step ri\ is expressed through u with the help of the continuity equa- 
tion: rii = —d(n u)/dx. To derive the energy conservation law we multiply Eq. (1) 
by du/dt. Then, using (||D, the first term in the right-hand side can be rewritten as 



2 dt 
n 2 h 2 d 



d(riQu) 
dx 
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_l_ d a V i u) a(n u) ^ ^ Q rema j n j n g terms can be combined into 



d(nQu) d(nQu) 



dx dt 



d^f^ d(n u) 2 \ ( d(n u) \2 
dxy b ° dx ) dx I 



, ; , , , . After bringing all these terms into the left-hand side 

and integrating over x from to oo (taking into account that u(0,t) = 0), Eq. (1) reduces 
to dE/dt = 0, where the energy E is given by 



E 



dx 



mn (x) / du\ 

2 vdty 



2 7T 2 ^ 2 , , 

i Vb + n (x / 



m 



d(iiou) 
dx 



(6) 



The expression (||) for the energy allows one to write down the corresponding Hamilto- 
nian. Treating the displacement u(x) as an operator u(x) and introducing the conjugate 
momentum operator p(x), [u(x),p(x')] = ih5(x — x')}, one gets 



n 



dx 



p 2 (x) 



2mno(x) 2 V m J \ dx 



(7) 



If the concentration is constant, the last term in the integrand takes the standard form! 
^mn Vp(^) 2 , where vp — iin h/m is the Fermi velocity. 



III. GENERAL FORM OF THE CORRECTION TO THE DENSITY OF STATES 

In this section we will derive a general expression for the correction, caused by the 
variation of the electron concentration n (x), to the tunnel density of states at the edge of 
the Luttinger liquid. The Hamiltonian (^) can be diagonalized by the transformation 

u(x) = Y,-r^=®A x )Q* ( 8 ) 
n (x) 



where ^(x) are the eigenfunctions of the equation 

2 



no{x) i 



'irh 



d 



(9) 



and f2 2 are the eigenvalues. We assume that $^ are normalized: / °° c?x$^(x)$ Ai '(x) = 
S^i, and turn to zero at the edge: 3^(0) = 0. As a result of the transformation (g) the 
Hamiltonian (|?p reduces to the system of harmonic oscillators 



P 2 mft 2 
2m ' 2 



(10) 



with frequencies Q^. Next, the transformation flD is applied to the operator which creates 
an electron at the edge x = 



*w P (-if ^) = «p(-E«A). (") 



where the coefficients a M are defined as 



(1//= [°° dx lM=. (12) 
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One can also express a M explicitly through (d^^,/dx)\ x= o- Dividing Eq. (|D by yrio(x) and 
integrating, one gets 



(13) 

x=0 

With Hamiltonian and having the form (|10D and ([TT|) , the calculation of the density of 
states z/(cj) reduces to the standard procedure! and results in 

v(u) = -Re / dte iuJt < tf(t)tft(0) >= -Re / dte iuJt e~ w{t) , (14) 

IV JO TV JO 

where the function W(t) is the sum over eigenmodes 

At this point we will make use of the assumption that the relative variation of the 
concentration rio(x) is small: no(x) = n + 5n(x); 5n <C n. For the average concentration, 
n, the eigenf unctions of Eq. (|9|), which satisfy the boundary condition, are oc sinkx, the 
corresponding eigenvalues being s 2 k 2 , where 



W(t)=Y,'^{l-e-^). (15) 



mvp 



is the plasmon velocity. Then the function W(t) takes the standard form 

W (t) = — r^(l- e - iskt )e- r ° k = —\J r -^^-), (17) 
Vf Jo k vf V r / 

where, as usually, the cutoff parameter e~ r ° h is introduced to insure the convergence of the 

integral at k — > oo. With Wo(£) given by fllTD , Eq. (|TJ|) reproduces the result (|1|). Treating 

the difference W(t) — Wo(t) as a perturbation and expanding the exponent in flTjj ) to the 

first order, one gets the following correction to the density of states 

6 



Su 



1 



-Re / dte wt l 
7r Jo \r + ist 



r 



s/v F 



W(t) - W (t) 



(18) 



This expression can be simplified if one uses the identity 

1 



(r + ist) s/vF 



dzz^ e- z{ro+ist) . 



(19) 



T(s/v F ) Jo 

Substituting flT9) into (fiW) and performing integration first over t and then over z, one gets 



5v 



r 



s/v F 



(20) 



T(s/v F ) V s 

The divergence of the integral at small k is canceled by the similar divergence in the first 
term. Further calculations require the knowledge of the eigenfrequencies For weak 
periodic modulation they are found in the next Section. 



IV. PLASMON SPECTRUM 



To find the eigenvalues of Eq. (|9|) it is convenient to introduce a function 



$ M = $ M \Jn (x)(n (x) + V m/7r 2 h 2 ). 



Then the equation for $ M takes the form 



d 2 <t>, 
dx 2 



1 / drio/dx 



n {x] 



4 \n (x) + Vqui / Ti 2 Ti 2 / (n (x) + Vom/iT 2 h 2 ) \\nhn (x) 



Vt^m \ 2 d 2 n /dx 2 



2no(x 



(21) 



(22) 



Let us specify the modulation profile: 5n(x) = — n M cos(ax + (p), where (p is the phase 
and a = 2-nja. For such a modulation and i%m ~ n the plasmon spectrum was studied 



numerically in Ref. [10]. If Um <C n the analytical solution can be obtained. After expanding 
the "potential energy" in (^) to the first order in um/u, the equation reduces to the 
conventional Mathieu equation^ 

d 2 $, 



dx 2 



+ 



— - cos(crx + lf) + \ 

Zi \ s 



(23) 



where the dimensionless parameter 6 is defined as 



_ 2n M 



n 



as] [s 2 / 2s 2 



(24) 



Since 9 is proportional to the relative modulation, we have 9 <C 1, so that the gap in the 
spectrum of eigenfrequencies is narrow. The center of the gap is determined by the 
condition fl^/s = a/2. This allows to simplify the expression for 9 by the substitution 
= os I '2. Then (^) takes the form 

n M VjA _ ( n M \ Vom/n 2 h 2 



fl = _f!ifl_I£: = f_f!i u 7 (25) 
2n V s 2 / V2nyn + \/ m/7r 2 ^ 2 ' V ; 

We emphasize again that 9 is proportional to the interaction strength Vo; this supports the 
obvious observation that there is no gap for non-interacting electrons. Note also that, in 
principle, there are higher-order gaps in the spectrum of Eq. (|23|) . They are centered at 
f2 M = psa/2 (p = 2,3...) with the widths proportional to (nj\//ra) p . However, they do not 



describe the higher-order gaps in the spectrum of Eq. (p2|), since (|23| ) was derived from ([22]) 
using the first-order expansion. 

As it was mentioned in the Introduction, the eigenmodes of Eq. (|23|), satisfying the 
condition $ M (0) = 0, can be of two types: localized and extended. 



A. Localized mode 



Substituting into (P3|) 

~ s (JX 

$ = 2n— V7e~ 7x sin — , (26) 
Up z 

where the prefactor insures the normalization, and neglecting the "non-resonant" terms 
oc sin(^ + (p), we get the following expressions for the frequency and the decay constant 
of a localized mode 

n= YV /l + 0cos ^~ y (l + ^cosJ (27) 
o9 

7 = sin 09. 
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We see that the localized mode exists only if sirup < 0, when 7 is positive. As if changes 
within the interval [n, 2n], the frequency Q moves from the lower to the upper boundary of 
the gap. The decay constant of the localized mode turns to zero as its frequency merges 
with the continuum. To calculate the contribution of the localized mode to Sv, we need the 
constant a. Substituting ( P7D into (|TBD, we obtain 

4 Fy 



a = — - 



a v n 



(28) 



B. Extended modes 



Let L be the normalization length. We search for solution of (^) in the form 



fll COs( + Kfj)x + Ip^j + (3 2 COS^ - K^X + 1p 2 



(29) 



The boundary condition requires that 



(5\ COS Ipl + /?2 COS ^2 = 0. 



(30) 



Substituting ( p9|) into fl23"D , we establish the following relation between phases 



(31) 



The system of coupled equations for coefficients /3 2 takes the form 



a 

2 ^ 



2n 

Pi + —P2 = 

2n 

+ —Pi = 0. 



(32) 



Solution of the system yields the plasmon spectrum 



A 2 

% = + AWk 2 s 2 + — , 



(33) 



where u = sa/2, and the width of the gap, A, is given by 



A = — — = 6u} , 



(34) 



so that the smallness of 9 guarantees that A/c^o <C 1. Parameter A takes the values ±1 
depending on the sign of 



A = sign^n^). 



(35) 



The plasmon spectrum (|33"D is shown in Fig. 1. The phase of modulation ip does not affect 
the spectrum, but it reveals itself in the parameter a M , which for given by (|29|) , takes 
the form 



a. 




2 Pi COS Tpi + Pi COS '02 



4i« (A 2 + /? 2 2 ) 1/2 



(36) 



The last factor in (|36"D can be expressed through <p using the relations 

Pi COS tpi + P 2 COS 1p2 Pi - Pi 



and fl3lD 



{Pi + piy 12 {Pi + p\yi 2 {P\ + pi + Wip 2 cos ^)V2 • 

Expressing the ratio Pi/ Pi with the use of Eqs. (03, (PHI), we finally obtain 



(37) 



/?! COS -01 + #2 COS 1p2 

W+W 1 ' 2 



SK. 



r. r, A2 . 2 

s K « + T Sln V 9 



1 + A COS <£> 



A/2 



1/2 



(38) 



Discreteness of is related to the final normalization length. Summation over /i reduces 
to integration with the help of the relation Sk^ = tt/L. 



V. ASYMPTOTES AND NUMERICAL RESULTS 



The resulting expression for 8v emergies after substitution of (|27|), (f28|), (|33|) and (|36|) 



into Eq. (|20|). It is convenient to present this expression by introducing the relative deviation 
of the frequency uo from the center of the gap 

Q = 2( " A " o) . (39) 

Then we obtain 

1 /r 6»y/^ / s \ dr-i 



s ^ = rW^){-f) W*' (40) 



where the dimensionless function F{Co) has the form 
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F{u) = 

and the functions fi and fi are defined as 



/- \ — — 1 1 £ cos 
7rsm<p(w - COS<p)«F --/i(^)H —J2{u) 



(41) 



/i(tD) = / tfo 



OO 



OC 



(d) — uVf 1 5 — - — 2 — (a) - X v Vv 2 + 1W 1 

ir + sin ip 



(42) 



/ 2 (a>) = dv X f (Q - KV^Tl) '* 1 . (43) 

J-oo [v z + sin <p) yv z + 1 v 7 

Similar to (p5|) , in ([421), (f!3|) X v = sign(v); the form (■■■) VF is defined only when the 
argument is positive and it is zero otherwise. 

We remind that the first term in F(u) is present only if siny? < 0. This term desribes 
a power-law anomaly in the density of states originating from the localized plasmon. The 
second and the third terms represent the contribution from the continuous plasmon spectrum 
outside the gap. 

In the integral fx two terms of the integrand almost cancel each other at large negative 
v, indicating that the modulation-induced correction to the spectrum is important only 
in the region of the gap. One can see that the integral fi divergies as v — > — oo, since 
for large negative v the integrand behaves as —(—v) v f . The origin of the divergency is 
that in (f42"D , (|43|) we have used the expansion of the spectrum in the vicinity of the gap: 



|Jfe-f| < or, so the integration should be cut at v ~ — 1/9. Note however, that the divergent 
part is frequency-independent and causes only a small correction to the background density 
of states z/(u;o). To study the frequency dependence one should subtract this divergent part 
by considering the difference fi = fi{oj) — /2(0), which converges. 

For large positive u (outside the gap) all three contributions to F(uj) behave as uj v f 
with prefactors depending on the phase of modulation ip. For integrals fi and fi the following 
asymptotic behavor in the limit £0 ^> 1 can be derived 

fi ~ ir\ sin tp\u v F , 

h = M&) - /a(0) « ^W^" 1 - ( 44 ) 
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In the opposite limit uj < 0, \u\ ^> 1 the contribution from the localized plasmon is absent 
and the integral f\ falls off as \uo\ v f . The only contribution to F(u) in this limit comes 
from f'2 which has the following asymptotics 

The frequency dependence of the correction 8v{£u), calculated numerically from Eq. (|4lD 
for various values of the phase if, is shown in Fig. 2. It can be seen that some of the curves 
have cusps at frequencies, corresponding to the boundaries of the gap u> = ±1. The cusps 
in the curves (p = ll7r/8 and <p = 37r/2 have their origin in the localized plasmon. 



VI. CONCLUSION 

The singularity in the tunnel density of states, caused by the periodic modulation of the 
electron concentration, should reveal itself as an anomaly in the differential resistance dl/dV 
at voltage eV = Huq = nhs/a, which depends on the period of modulation and the plasmon 
velocity. We considered the case of a weak sinusoidal modulation. If higher harmonics 
with periods a/p are present in n (x), the gaps in the plasmon spectrum at u — pu (and 
correspoding anomalies in dl/dV at eV = np%s/a) should emerge. 

The strength of the anomaly is governed by the dimensionless parameter 9 which 
we assumed to be small. To sense the magnitude of the anomaly it is convenient to relate 
Sis to the density of states at u = uq in the absence of modulation: 8v/v(ujq) = 9 S ' VF F{Co). 

In our calculations we assumed the interactions to be short-range. This implies that the 
interaction radius ~ r should be much smaller than the period of modulation a. 

As a final remark, we considered a single-channel Luttinger liquid. Increasing of the 
number of channels would result in the new anomalies in the tunnel density of states, caused 
by resonances within each channel. Besides, modulation would couple the plasmons from 
different channels, thus causing additional anomalies. 
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FIGURES 

FIG. 1. Plasmon spectrum of a periodically modulated Luttinger liquid. 

FIG. 2. The normalized correction to the density of states, 5v /v(ujq)6 s I Vf , is plotted as a func- 
tion of the dimensionless frequency & for ^ = 1.2 and different values of the phase of modulation 
if. (a) p = (solid curve), p = it/ 4 (long-dashed curve), if = tt/2 (dash-dotted curve), (b) 
(p = 37r/4 (solid curve), ip = it (long-dashed curve), ip = ll7r/8 (dash-dotted curve), (c) ip = 3ir/2 
(solid curve), ip = 77r/4 (long-dashed curve), ip = 157r/8 (dash-dotted curve). 
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